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In a recent paper it was shown that if a Hamiltonian H has an unbroken VT symmetry, then it 
also possesses a hidden symmetry represented by the linear operator C. The operator C commutes 
with both H and VT. The inner product with respect to CVT is associated with a positive norm 
and the quantum theory built on the associated Hilbert space is unitary. In this paper it is shown 
how to construct the operator C for the non-Hermitian PT-symmetric Hamiltonian H — |p 2 + 
|a; 2 + íex 3 using perturbative techniques. It is also shown how to construct the operator C for 
H — |p 2 + t;X 2 — ex 4 using nonperturbative methods. 
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O ■ I. INTRODUCTION AND BACKGROUND 

Z: 

' It was observed in 1998 [1] that with properly defined boundary conditions the Sturm-Liouville difïerential equation 

ÇSJ , eigenvalue problem associated with the non-Hermitian Hamiltonian 

i— T H=p 2 + x 2 (ix) u O>0) (1.1) 

, exhibits a spectrum that is real and positive. It was argued in Ref. [1] that the reality of the spectrum of H is a 
' consequence of the unbroken VT symmetry of H . A complete proof that the spectrum of H is real and positive was 
given by Dorey et al [2] . 

By VT symmetry we mean the following: The linear parity operator V performs spatial reflection and thus reverses 
çvq , the sign of the momentum and position operators: VpV^ 1 = —p and VxV^ 1 = —x. The antilinear time-reversal 
operator T reverses the sign of the momentum operator and performs complex conjugation: TpT~ L = —p, TxT^ 1 = 
x, and TiT^ 1 = —i. The Heisenberg àlgebra, [x,p] = i, which is fundamental in quantum theory because it embodies 
O ,' the uncertainty principle, is invariant under the action of the operators V and T separately. The non-Hermitian 
Hamiltonian H in (1.1) is not symmetric under V or T separately, but it is invariant under their combined operation; 
' such Hamiltonians are said to possess space-time reflection symmetry (VT symmetry). 

We say that the VT symmetry of a Hamiltonian H is not spontaneously broken if the eigenfunctions of H are 
simultaneously eigenfunctions of the VT operator. It is difficult to prové that the VT symmetry of a given Hamiltonian 
is not spontaneously broken, but if this is the case, then it is easy to show that the spectrum is entirely real [3]. 

Space-time reflection (VT) symmetry is a weaker condition than Hermiticity in the following sense. For many 
different Hermitian Hamiltonians, such as H = p 2 + x 4 , H = p 2 + x 6 , H = p 2 + x s , and so on, we can construct 
1 infinite classes of non-Hermitian PT-symmetric Hamiltonians H = p 2 + x 4 (ix) v ' , H = p 2 + x & (ix) u , H = p 2 + x 8 (ix) I/ , 
and so on. So long as the parameter v is real and positive (v > 0), the VT symmetry of each of these Hamiltonians 
is not spontaneously broken and the spectrum is entirely real [3] . 

Showing that the Sturm-Liouville problem associated with a non-Hermitian 'PT-symmetric Hamiltonian has a 
positive real spectrum is mathematically significant, but it does not have any obvious relevance to physics. To show 
that a Hamiltonian can serve as the basis for a theory of quantum mechanics it is necessary to demonstrate that the 
Hamiltonian acts on a Hilbert space that is endowed with an inner product whose associated norm is positive definite. 
Only then can one say that the theory has a probabilistic interpretation. Furthermore, it must be shown that the 
theory is unitary (probability must be conserved in time). Since the publication of Ref. [1] it has been believed that 
the Hamiltonians in (1.1) could not be the basis for a physical theory because they are non-Hermitian. Indeed, the VT 
norm is not positive definite and this appears to present interpretational problems in developing a quantum theory 
based on PT-symmetric Hamiltonians. Many papers have been published that discuss this apparent shorteoming of 
non-Hermitian PT-symmetric Hamiltonians [4]. 

In a recent letter it was shown how to overcome these problems [5] . This letter demonstrates that any Hamiltonian 
that possesses an unbroken VT symmetry also has a hidden symmetry. This new symmetry is represented by the 
linear operator C, which commutes with both the Hamiltonian H and the VT operator. In terms of C one can 
construct an inner product whose associated norm is positive definite. Observables exhibit CVT symmetry and the 
dynamics is governed by unitary time evolution. Thus, PT-symmetric Hamiltonians give rise to new classes of fully 
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consistent complex quantum thcorics. Thesc ncw quantum theories are extensions of conventional Hcrmitian quantum 
mechanics into the complex domain. The novelty of these theories is that the inner product is not specified prior to 
and independently of the Hamiltonian. Rather, the inner product is determined by the Hamiltonian itself. Thus, in 
such theories the norm and hence the notion of probability is dynamically incorporated in the Hamiltonian. 

The purpose of the present paper is to present an explicit calculation of C for two nontrivial Hamiltonians. First, 
we consider the case of the PT-symmetric Hamiltonian 

H = ^p 2 + ^x 2 +iex 3 , (1.2) 

for which we give a perturbative calculation of the operator C correct to third order in powers of e. Second, we 
calculate C for the Hamiltonian 

H=\p 2 + \J ex\ (1.3) 

for which ordinary perturbative methods arc ineffcctive and nonperturbativc methods must be used. The organization 
of this paper is straightforward. In Sec. II we review the formal construction, first presented in Ref. [5], of the C 
operator. In Sec. III we calculate C for the Hamiltonian in (1.2) and in Sec. IV we calculate C for the Hamiltonian in 
(1.3). 

II. FORMAL DERIVATION OF THE C OPERATOR 

In this section we present a formal discussion of PT-symmetric Hamiltonians and we show how to construct the C 
operator. In general, for any PT-symmetric Hamiltonian H we must begin by solving the Sturm-Liouville differential 
equation eigenvalue problem associated with H: 

H<t> n {x) = E n 4> n (x) (n = 0, 1, 2, 3, ■■■)■ (2.1) 

For Hamiltonians like those in (1.1 - 1.3) the differential equation (2.1) must be imposed on an infmite contour 
in the complex- x plane. For large |x| the contour lies in wedges that are placed symmetrically with respect to the 
imaginary-x axis. These wedges are described in Ref. [1]. The boundary conditions on the eigenfunctions are that 
4>(x) — > exponentially rapidly as \x\ -^coon the contour. For H in (1.2) the contour may be taken to be the real-x 
axis, but for H in (1.3) the contour lies in the two wedges — 7r/3 < argrr < and — 7r < argx < — 27r/3. It is not 
possible to solve the differential equation (2.1) analytically for the two Hamiltonians (1.2) and (1.3) considered in this 
paper but we have solved it numerically to very high aceuracy for the first ten eigenfunctions and eigenvalues. As 
mentioned above, the eigenvalues are all real and positive and are nondegenerate. 

For all n, the eigenfunctions <j> n (x) are simultaneously eigenstates of the VT operator: VT(j> n (x) — X n <j> n (x). 
Moreover, because (VT) 2 — 1 and VT involves complex conjugation, it follows that |A„| = 1. Thus, A„ = e M " is a 
pure phase. For each n this phase can be absorbed into </>„ by the multiplicative rescaling <j> n — > e~ M ™/ 2 0„, so that 
the new eigenvalue of VT is unity: 

VT<j> n (x)=<j> n (x). (2.2) 

Next, we observe that there is an inner product, called the VT inner product, with respect to which the eigenfunc- 
tions 4> n (x) for two different vàlues of n are orthogonal. For the two functions f(x) and g(x) the VT inner product 
(/, g) is defined by 

(f,g)= í dx [VTf(x)}g(x), (2.3) 
Jc 

where VTf(x) = [/(— x)]* and the contour C lies in the wedges described above. For this inner product the associated 
norm (/, /) is independent of the overall phase of f(x) and is conserved in time. (Phase independence is required 
because ultimately we must construct a space of rays to represent quantum mechanical states.) The proof that 
eigenfunctions <f> n (x) corresponding to different vàlues of n are orthogonal with respect to this inner product is trivial 
and follows directly from the differential equation (2.1) using integration by parts. 

We then normalize the eigenfunctions so that \(4> ni <j) n )\ = 1 and we discover the apparent problem with using a 
non-Hermitian Hamiltonian. While the eigenfunctions are orthogonal, the VT norm is not positive definite: 

{4>mAn) = {-l) n 5 m ,n (m, n = 0, 1, 2, 3, ■ ■ ■ )■ (2.4) 
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Despite the fact that this norm is not positive definite, the eigenfunctions are complete. For real x and y the 
statement of completeness in coordinate space is 1 

52(-ï) n Mx)<t>n(v) = 8(x-y). (2.5) 

n 

This is a nontrivial result that has been verified numcrically to extremely high accuracy [6]. Using (2.4) we can vcrify 
that the sum in (2.5) is the position-space representation of the unity operator: 



dy S(x — y)S(y — z) = 5(x — z). (2.6) 
We can also express the Hamiltonian H and the Green's function G(x, y) in the coordinate space representation: 
H{x,y) = Y,{-^ n En4>n{x)My) and G(x, y) = ^WÍ-M^Mv)- ( 2 -7) 

n n 

The Green's function G(x, y) satisfies the inhomogeneous differential equation 

HG(x,y) = S(x-y), (2.8) 

which states that the Green's function is the inverse of the Hamiltonian operator. 

In addition, we can construct the parity operator V in terms of the energy eigenstates. In position space 

V(x,y) = S(x + y) = ^(-l)> n (x)^ n (-y). (2.9) 

n 

Again, using (2.4) we can see that the square of the parity operator is unity. 

Finally, we construct the linear operator C that expresses the hidden symmetry of the Hamiltonian H . The position- 
space representation of C is 

C(x,y) =Y,<t>n{x)<t>n{y). (2.10) 

n 

The properties of the operator C are easy to verify using (2.4). First, like the parity operator, the square of C is unity: 

dyC(x,y)C(y,z) = d(x-z). (2.11) 



Second, the eigenfunctions </>„ (x) of the Hamiltonian H are also eigenfunctions of C and the corresponding eigenvalues 
are (-1)": 



/ 



dyC{x,y)^ n {y) = {-l) n cj> n {x). (2.12) 



Third, the operator C commutes with both the Hamiltonian H and the operator VT. Note that while the opcrators 
V and C are unequal (the parity operator V is real, while the operator C is complex), both V and C are square roots 
of the unity operator 5{x — y). Last, the operators V and C do not commute. Indeed, CV = (VC)*. 

The operator C docs not exist as a distinct entity in conventional Hermitian quantum mechanics. Indeed, we will 
see that as the parameter e in (1.2) and (1.3) tends to zero the operator C becomes idcntical to V . Thus, in this limit 
the CVT operator becomes T. This verifies that for symmetric Hamiltonians in Standard quantum mechanics CVT 



l· It is important to remark here that the argument of the Dirac delta function in (2.5) must be real because the delta function 
is only defined for real argument. This may seem to be in conflict with the earlier remark in this section that the Schrodinger 
equation (2.1) must be solved along a contour that lies in wedges in the complex- x plane. To resolve this apparent conflict we 
specify the contour as follows. We demand that the contour lie on the real axis until it passes the points x and y. Only then 
may it veer off into the complex-^ plane and enter the wedges. This choice of contour is allowed because the wedge conditions 
are asymptotic conditions. The positions of the wedges are determined by the boundary conditions. 
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symmetry and Hermiticity coincide and CVT invariance can be viewed as the natural complex extension of the usual 
Hermiticity condition. 

We can now define an inner product (f\g) whose associated norm is positive: 



(f\g) 



= I dx[CVTf{x)]g{x). (2.13) 



The CVT norm associated with this inner product is positive because C contributes —1 when it acts on states with 
negative VT norm. To verify that this norm is positive definite we expand an arbitrary function f(x) as a linear 
combination of cigenfunctions of the Hamiltonian H: 



Then, the CVT norm of f(x) is 

/■oo oo 
dx[CVTf(x)]f(x) = J2\cn\ 2 , 



which is positive unless f(x) = 0. The CVT norm is time independent because the CVT operator commutes with the 
Hamiltonian H and thus the theory is unitary. Using the CVT conjugate, the completeness condition (2.5) becomes 

J2iCVTMx)}My) = S(x - y). (2.14) 



III. PERTURBATIVE CALCULATION OF C IN A PT-SYMMETRIC CÚBIC THEORY 

In this section we use perturbative methods to calculate the operator C(x, y) for the Hamiltonian H = ip 2 + ^x 2 + 
iex 3 . We perform the calculations to third order in perturbation theory. We begin by solving the Schròdinger equation 

-\t>n{x) + ^x 2 (f) n (x) + iex 3 (t) n (x) = E n (j) n {x) (3.1) 

as a series in powers of e. 

The perturbative solution to this equation has the form 

4> n {x) = ^ i/4 ^^_ e-^ 2 [H n (x) - iP n {x)e - Q n (x)e 2 + iR n (x)e 3 ] , (3.2) 

where H n (x) is the nth Hermite polynomial and P n (x), Q n (x), and R n (x) are polynomials in x of degree n + 3, n + 6, 
and n + 9, respectively. These polynomials can be expressed as linear combinations of Hermite polynomials: 



4 



13 3 1 

p n{x) = —H n+3 (x) + -(n + í)H n+1 (x) - -r^Hn-xix) - -n(n - l)(n - 2)i?„_ 3 (a;) 



J_ií„ +6 (x) + -^(4n + 7)H n+i {x) + ^(7n 2 + 33n + 27)H n+2 (x) 



1 



1 



+ -n(n - l)(7n 2 - 19n + l)ií„_ 2 (a;) + -n(n - l)(n - 2)(n - 3)(4n - 3)i?„_ 4 (a;) 

o 5 

+ - l)( n - 2)(n - 3)(n - 4)(n - 5)JT n _ 6 (a;), 



82944 



#«+9(2) 



3072 



(2n + 5)iï„ +7 (x) 



7680 



(80n 2 + 465n + 549)7í„ + 5(x) 



+ 



1 



6912 
3 



(488n 3 + 3639n 2 + 9832n + 7506)H n+3 (x) 



+ — (20n 4 - n 3 + 203n 2 + 408n + 228)H n+1 (x) 

--^n(20n 4 + 81n 3 + 326n 2 + 81n + 44)ií„_ 1 (a;) 
64 

- l)(n - 2)(488n 3 - 2175n 2 + 4018n - 825)iJ„_ 3 (:r) 

864 

~2Í0 n( - n - l)(n - 2)(n - 3)(n - 4)(80n 2 - 305n + 164)iJ n _ 5 (a;) 
~^n{n - l)(n - 2)(n - 3)(n - 4)(n - 5)(n - 6)(2n - 3)ií n _ 7 (x) 
-^2 n (" - - 2 )( n - 3 )( n - 4 X™ - 5 )( n - 6)(n - 7)(n - 8)H n _ 9 (x). 



The energy i? n to ordcr e 3 is 



E n =n+i + i(30n 2 + 30n+ll)e 2 + 0(e 4 ). 
2 o 



The expression for (f> n (x) must be 'PT-normalized according to (2.4) so that its square integral is (—1)™: 



f 

J — ( 



dx [0„(x)] 2 = (-l)" + O(e 4 ). 



This determines the value of a n in (3.2): 



a n = 1 



144 



(2n + l)(82n 2 + 82n + 87)e 2 + 0(e 4 ). 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



We calculate the operator C(x,y) = Y^=o ( t ) n{x)4> n {y), which is given formally in (2.10), by directly substituting 
the wave functions (j) n {x) in (3.2). We then use the completeness relation for Hermite polynomials, 



1 00 1 

^ {X2+V2) E ^ H ^) H n{y) = S(x - y), 



n=0 



2 n nl 



(3.7) 



to evaluate the sum. We also need to use the following identities satisfied by the Hermite polynomials: 



xH n (x) = -H n+1 (x) +nH n _ 1 (x), 

H';{x) = 2xH' n {x)-2nH n {x), 
H' n {x) = 2nH n _ 1 (x). 



To third order in e the result is 
C(x,y) 

^J^ + l3^- 48 ^J^ + (- 8 ^ + 64 ^ 



, . /4 9 3 „ d 
l - le \zdx-z +2xy dx- 



8 d e 8 d 4 fn 2 2 d 2 
9^ + 3^ + (2 ^- 12 te 



8 2 2 176\ d 5 /4 3 3 , \ d 3 ., . 

3 y 



<9 



(3.8) 

32 d 9 16 a 7 
8Ïdx^ + Y xy dxï 

+ 0(e A )\s(x + y). (3.9) 



Hcnce, the coordinate-space representation of the operator C{x 1 y) is expressed as a derivativa of a Dirac delta function. 
From this expression for C(x, y) we can verify the following properties: First, to order e 3 the operator C{x, y) satisfics 
(2.11). That is, 



f 

J — ( 



dy C(x, y)C(y, z) = S(x - z) + 0(e 4 ). (3.10) 



Second, to order e 3 the operator C(x,y) satisfics (2.12); the wave functions <j) n {x) are eigenstates of C(x,y) with 
eigenvalue (—!)"• That is, 



f 

J —i 



dyC(x,y)My) = (-l) n Mx) + 0(e 4 ). (3.11) 



Third, in the limit as e — > 0, the operator C(x,y) becomes the coordinate-space representation of the parity operator 
V(x,y) = S(x + y). 

There is a somewhat simpler way to express the operator C(x, y). The derivative operator in (3.9) that is acting on 
S(x + y) can be exponentiated so that to order e 4 (and not just e 3 ) we have 

C(x, y) = e- UA -^ B 5{x + y) + 0(e 5 ), (3.12) 
where the derivative operators A and B are given by 

A--—-2x—x 

3 dx 3 dx 

„ 128 d 5 40 d 3 n 2 d 2 „„ d ,„ ,„, 

B = ¥^-ïV + & V- 3 ^ (3 ' 13) 

We have applied the procedure used above to calculate C(x,y) to cvaluate the parity operator V(x,y). That is, 
we have substituted the eigenfunctions <j> n (x) in (3.2) into the formal sum in (2.9). We find that to each order in 
powers of e the summation vanishes except for the leading term (the coefficient of e°). Thus, we obtain the result that 
V(x, y) = S(x + y) + 0(e 4 ). This is not a new result, but it provides a useful check of the aceuracy of our calculations. 
Similarly, we have evaluated the sum in (2.5) and we obtain the trivial result S(x — y) + 0(e 4 ). We have also evaluated 
the expression in (2.7) for the Hamiltonian in coordinate space and we find (as expected) that the coefficient of e k in 
the summation vanishes for k > 1 and we get 

H{X ' V) = {~l& + + ^ 5{X - y) + ° (e4) - 

We have again applied the procedure for calculating C(x,y) to evaluate the Green's function G(x,y) in (2.7). Sub- 
stituting the eigenfunctions <j) n {x) in (3.2) into (2.7) and performing the summation gives the perturbative expansion 
of the Green's function: 

G(x, y) - G (x, y) - iG x {x, y)e - G 2 (x, y)e 2 + iG 3 (x, y)e 3 + 0(e 4 ). (3.14) 
The zeroth-order Green's function satisfies the differential equation 

\l^ + \ x2 ) G ^ y) = 5{x - y) - (3 ' 15) 
The solution to this equation is 

G (x,y) - 6(x - y)D_ 1/2 (xV2)D_ 1/2 (-yV2) + 9(y - x)D_ 1/2 (-xV2)D_ 1/2 (yV2), (3.16) 
where D v (x) is the parabòlic cylindcr function and 9(x) is the step function defined by 

( {x < 0), 

9(x) = \ \ (x = 0), (3.17) 
[ 1 (x> 0). 

Note that Go(x,y) is a symmetric function of x and y. 
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The first-order contribution to the Green's function satisfies the differential equation 

1 8 2 1 



. );ií2 +-x 2 )G 1 (x 1 y)=x 3 G Q (x,y) (3.18) 



and the solution to this equation is 



Gi(x, y) = -Ux 2 ^--x + y 2 ^- - y ) G (x,y). (3.19) 



3 \ dx dy 
The second-order contribution to the Green's function satisfies 

1 d 2 1 



. 2i)y , 2 x 2 )G 2 (x 1 y)=x 3 G 1 (x,y) (3.20) 

and the solution to this equation is 

1 / 8 8 \ 2 7 r°° 

G 2 (x,y) = -^ 2 — -x + y 2 — -yj G Q (x,y) + - J dz z 4 G (z,x)G Q (z,y). (3.21) 

The third-order contribution to the Green's function satisfies 

^1 8 2 1 

and the solution to this equation is 



^ + ^x 2 )G 3 (x,y)=x 3 G 2 (x,y) (3.22) 



G 3 (x,y) = ~ 



( 5 o 1 o „ 56 4 112\ d 25 7 1 c 112 , 

— x s H x 2 y b H x 4 H — H x 7 x b y x A 

V36 12 y 15 5 J dx 36 12 y 15 

Go(x,y) 



1 K o 56 „ 112\ 9 25 - 1 í? 112 



36 r + Ï2 x6y2 + Ï5 y4 + —)d-y + 36 y7 - 12^ " 



7/8 8 \ í°° 

--^c 2 —-x + y 2 —-y)J dzz 4 G (z,x)G (z,y). (3.23) 



IV. NONPERTURBATIVE CALCULATION OF C IN A PT-SYMMETRIC QUARTIC THEORY 

In this section we explain briefly the nonpcrturbative methods that must be used to calculate the operator C(x,y) 
for the Hamiltonian H = ^p 2 + \x 2 — ex 4 . We follow the approach taken in Ref. [7] , in which nonperturbative methods 
were used to calculate the one-point Green's function for this Hamiltonian. 

A. Failure of Perturbation Theory 

We begin by explaining why perturbation theory fails to produce the operator C(x,y). Following the approach 
taken in Sec. III, we expand the solution to the Schrodinger equation 

-\<i>n{x) + \x 2 c^ n {x) - ex 4 cj) n (x) = E n cj) n (x) (4.1) 

as a series in powers of e: 

M*) = + We] + 0(e 2 ), (4.2) 

where H n (x) is the nth Hermite polynomial and P n (x) is a polynomial in x of degree n + 4. The polynomial P n {x) 
is a linear combination of Hermite polynomials: 

P n {x) = ^H n+4 {x) + l(2n + 3)H n+2 (x) - l -n(n - l)(2n - l)H n _ 2 (x) - -n(n - l)(n - 2)(n - 3)H n _ 4 {x). (4.3) 
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The energy E n to order e is 

E n = n+^-'^(2n 2 + 2n + l) e + 0(e 2 ). (4.4) 
We must also VT normalize the expression for <p n (x) according to (2.4) so that its square integral is (—1)™: 

/oc 

dx [0„(x)] 2 = (-l) n + O(e 2 ). (4.5) 
-oo 

This determines the value of a n in (4.2). The result is very simple; to order e we have 

a„ = l + 0(e 2 ). (4.6) 

Finally, we substitute 4> n (x) in (4.2) into (2.10) and use the identity in (3.7). However, we obtain the trivial result 
that only the leading term (zeroth-order in powers of e) survives. More generally, we can show by a parity argument 
that the coefficicnts of all highcr powers of e vanish. Thus, we get the (wrong) result that 

C(x,y) = S(x + y) (WRONG!). (4.7) 

We know that this result is wrong because the operator C{x 1 y) is complex and the result in (4.7) is real. An alternativc 
way to see this is to note (4.7) imples that C(x, y) and V(x, y) coincide; but in this PT-symmetric theory, C(x, y) and 
V(x,y) are distinct operators. We will see that the difference between C(x,y) and V(x,y) is a nonperturbative term 
of order e _1 /( 3£ \ which is smaller than any integer power of e. 



B. Nonperturbative Analysis 

We will now show how to perform a nonperturbative analysis of the Schròdinger equation (4.1). We decompose the 
eigenfunction (f> n (x) into its perturbative part on the right side of (4.2) and a nonperturbative part: 

<j> n (x) = <j>%^(x)+cf>™^(x). (4.8) 

The nonperturbative part of <f) n (x) is cxponcntially small compared with the perturbative part, but these two contri- 
butions can be easily distinguished because for real argument x, one is real while the other is imaginary. 

Following the WKB analysis in Ref. [7], we break the real-x axis into three regions: In region I, where |x| <C e -1 / 4 , 
we have 

i n 

cpr\x) ~ ^^D„(x^), 

7T ' V 7"ï<' 

„o„pcrt (x) _ lbnCn ( x ^2), (4.9) 

where the coemeient of D n is taken from (4.2) and the coeficient ib n of C n will be determined by asymptotic 
matching. Note that for nonnegative integer index the parabòlic cylinder function D n is expressed in terms of a 
Hermitc polynomial H n as 

D n (xy/2) = 2-"/ 2 e-^X(z). (4.10) 

Also, for nonnegative integer index the functions D n and C n are a pair of linearly independent solutions to the 
parabòlic cylinder equation. They can be expressed in terms of parabòlic cylinder functions as follows: 

D n (z) = -= [i n D- n -!{iz) + (-i) n D- n -!(-iz)] , 

X/ZTT 



C n (z) = -1= [i n D_„_i(w) - (-i) n D_ n _i(-«)] . (4.11) 

VZTT 



In region II, where 1 -C \x\ <C e -1 ^ 2 , we can obtain the eigenfunction using WKB theory. We write the Schròdinger 
equation (4.1) in the form <j)' r [(x) — u) n (x)(p n (x) where, to leading order in e, we have u n {x) — —2ex 4 + x 2 — 2n — 1. 
Then, for positive x the physical-optics WKB approximation reads 
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^ ert (x)~/„[o; n (x)]- 1 / 4 exp 

C o„pcrt (x) _ 5n[[Jíi(x)r l/4 exp 



ÚS ^üJ n (s) 

/ ds y/u n (s) 
J xi 



(4.12) 



where the constants /„ and g n will be determined by asymptotic matching. The lower endpoint of integration, 
X\ = \/2n + 1, is the approximatc location of the inner turning point. 

In rcgion III x is near the outer turning points at ±l/-\/2e. For positive x we define the variable r by x = 
x 2 (l — 2 1 / 3 e 2 / 3 r) , where x 2 = 1/V2e. The condition that x is near x 2 is that r <C e~ 2 / 3 . In this region the 
Schròdingcr cquation becomes an Airy equation in the variable r: <^(r) = r<p n {r). The solution in this region reads 



0P crt (O ^ fcnBi(r), 

nonpert (r) „ _ ifen Ai(r), 



(4.13) 



where Ai(r) and Bi(r) are the exponcntially decaying and growing Airy functions for large positive r. The fact that 
the same coefhcient h n multiplies both Bi and Ai is a nontrivial result that is established in Rcf. [7] . 

By asymptotically matching the solutions in regions I and II and the solutions in regions II and III we obtain the 
formula for the coefHcient of the nonperturbative part of the solution in (4.9): 



K = - 



V2d 



(4/e) 



1+1/2 



e 3 £ . 



(4.14) 



Finally, using the wave function in region I we can construct the operator C(x,y) according to (2.10): 

oo 

C{x,y) = ^<p n {x)<j) n {y) 

oo 

= E [€r\^r\y) + € CIt (xWn onpcit (y) + C onpcrt (^)^ ort (y) + r n ODpcit ^W n onpcit (y)] ■ (4.15) 



n=0 



The first sum in this equation gives S(x + y) to all orders in powers of e as explained above in Subsec. IVA. The last 
sum is negligible compared with the second and third sums. We thus obtain 



C(x, y) - 6(x + y)- i^Zfe e~h ^ ^(-4/e) n [D n (xV2)C n (yV2) + C n {xs/2)D n {ys/2) 



n=0 



(4.16) 



where C n and D n are defined in (4.11). Observe that the correction to the delta function (that is, the difference 
bctween the V operator and the C operator) is nonperturbative; it is exponcntially small and imaginary. 
The summation in (4.16) can be converted to a doublc integral: 



C(x,y) =d(x + y) + i 




2 e -è e è(- 2 +^ 2 ) 



d_ 

dx 



Jo Jo 



ds 



ViTs- 



: exp 



(2y/2s/e cos 6 — ix — isy 
1 + s 2 



+ (x <-» y) 



(4.17) 



This is the lcading-order nonperturbative approximation to the coordinate-space representation of the operator C. 
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